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The thermoelectric response of SND configuration is considered within the generahzed Ginzburg- 
Landau theory for a homogeneous admixture of s-wave and d-wave superconductors. The resulting 
thermopower AS{T, 6, tc) — Sp{9, tc) — B{9, tc){Tcs — T) is found to strongly depend on the relative 
phase 9 = (f)s — <j)d between the two superconductors {tc = Tcd/Tcs where Tcs and T^d are the 
corresponding critical temperatures). Two independent mechanisms are shown to contribute to the 
peak value Sp{0,tc). One, based on the charge imbalance between the quasiparticles and Cooper 
pairs (described by the corresponding chemical potentials, /i, and fip) due to the normal metal 
insert, results in a pronounced maximum of the peak near 6 = ±7r/2 (the so-called s ± id mixed 
pairing state) for two identical superconductors with Ted = Tcs- This mechanism can be realized 
in a d-wave orthorhombic sample (like YBCO) with twin boundaries which are represented by 
tetragonal regions of variable width, with a reduced chemical potential. Another mechanism (not 
related to the charge imbalance effects) occurs when two different superconductors with T^d 7^ Tcs 
are used for SND junction. It gives rise to Sp{0,tc) ex. 1 — tc and can be realized via the junction 
comprising an s-wave low- Tc superconductor (like Pb) and a d-wave high-Tc superconductor (like 
orthorhombic YBCO). The experimental conditions (based on the previous experience with SNS 
junctions) under which the predicted behavior of the induced differential thermopower can be mea- 
sured are discussed. 
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I. INTRODUCTION 

During the last few years the order parameter sym- 
metry has been one of the intensively debated issues in 
the field of high-Tc superconductivity (HTS). A number 
of experiments points to its d^2_y2-wa,Ye character [Q. 
Such an unconventional symmetry of the order param- 
eter has also important implications for the Josephson 
physics because for a d-wave superconductor the Joseph- 
son coupling is subject to an additional phase depen- 
dence caused by the internal phase structure of the wave 
function. The phase properties of the Josephson effect 
have been discussed within the framework of the gener- 
alized Ginzburg-Landau (GL) |^ as well as the tunnel- 
ing Hamiltonian approach |3j. It was found Q that the 
current-phase relationship depends on the mutual orien- 
tation of the two coupled superconductors and their in- 
terface. This property is the basis of all the phase sensi- 
tive experiments probing the order parameter symmetry. 
In particular, it is possible to create multiply connected 
d-wave superconductors which generate half-integer flux 
quanta as observed in experiments j^. Various inter- 
esting phenomena occur in interfaces of d-wave super- 
conductors. For example, for an interface to a normal 
metal a bound state appears at zero energy giving rise to 



a zero-bias anomaly in the /-^-characteristics of quasi- 
particle tunneling while in such an interface to an 
s-wave superconductor the energy minimum corresponds 
to a Josephson phase different from or tt. By sym- 
metry, a small s-wave component always coexists with a 
predominantly d-wave order parameter in an orthorhom- 
bic superconductor such as YBCO, and changes its sign 
across a twin boundary ||^. Besides, the s-wave and d- 
wave order parameters can form a complex combination, 
the so-called s ± id-state which is characterized by a local 
breakdown of time reversal symmetry T either near sur- 
faces or near the twin boundaries represented by 
tetragonal regions with a reduced chemical potential p3[ . 
Both scenarios lead to a phase difference of ±7r/2, which 
corresponds to two degenerate states Moreover, 
the relative phase oscillations between two condensates 
with different order parameter symmetries could mani- 
fest themselves through the specific collective excitations 
("phasons") 

At the same time, a rather sensitive differential tech- 
nique to probe sample inhomogeneity for temperatures 
just below Tc, where phase slippage events play an im- 
portant role in transport characteristics has been pro- 
posed and successfully applied for detecting 
small changes in thermopower of a specimen due to the 
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deliberate insertion of a macroscopic S* TVS' junction made 
of a normal-metal layer N, used to force pair breaking 
of the superconducting component when it flows down 
the temperature gradient. In particular, a carrier-type- 
dependent thermoelectric response of such a SNS con- 
figuration in a C-shaped BixPbi-xSr2CaCu20y sample 
has been registered and its A-shaped temperature behav- 
ior around Tc has been explained within the framework 
of GL theory @,[9|. 

In the present paper, we consider theoretically the case 
of SND junction and discuss its possible implications for 
the above-mentioned type of experiments. The paper is 
organized as follows. In Section II we briefly review the 
experimental results for SNS configuration (with both 
holelike and electronlike carriers of the normal-metal N 
insert) and present a theoretical interpretation of these 
results, based on GL free energy functional. The cru- 
cial role of the difference between the quasiparticle /ig 
and pair /ip chemical potentials in understanding the ob- 
served phenomena is emphasized. In Section III, extend- 



ing the early suggested |10[|13| GL theory of an admixture 
of s-wave and d-wave superconductors by taking into ac- 
count pair- breaking effects with ^ fj,p, we calculate 
the differential thermopower AS* of SND configuration 
near Tc- The main theoretical result of this Section is 
prediction of a rather specific dependence of AS* on rel- 
ative phase shift 9 = (ps — i'd between the two supercon- 
ductors. Two independent mechanisms contributing to 
the peak value of the thermopower are discussed. One, 
based on the charge imbalance between the quasiparti- 
cles and Cooper pairs due to the normal metal insert, is 
discussed in Section IIIA. It results in a pronounced max- 
imum of the peak near 9 ~ ±7r/2 (the so-called s ± id 
mixed pairing state) for two identical superconductors 
with Ted — Tcs = Tc- This mechanism can be realized, 
e.g., in a d-wave orthorhombic sample (like YBCO) with 
twin boundaries which are represented by tetragonal re- 
gions of variable width, with a reduced chemical poten- 
tial. Another mechanism (not related to the charge im- 
balance effects), discussed in Section IIIB, occurs when 
two different superconductors with T^d ^ T'cs are used 
for SND junction. This situation can be realized for 
an s-wave low- Tc superconductor (like Ph) and a d-wave 
high- Tc superconductor (like orthorhombic YBCO). 



II. SNS CONFIGURATION: A REVIEW 

A. Experimental setup and main results 

Before turning to the main subject of the present pa- 
per, let us briefly review the previous results concerning a 
carrier-type-dependent thermoelectric response of SNS 
configuration in a C-shaped BixPbi-xSr2CaCu20y 
sample (see Ref. Q for details). The sample geometry 
used is sketched in Fig.l, where the contact arrangement 
and the position of the sample with respect to the tem- 



perature gradient V^jT is shown as well. Two cuts are 
inserted at 90° to each other into a ring-shaped supercon- 
ducting sample. The first cut lies parallel to the applied 
temperature gradient serving to define a vertical symme- 
try axis. The second cut lies in the middle of the right 
wing, normal to the symmetry axis, separating an s-wave 
superconductor (S" = 5*) from another s-wave supercon- 
ductor {S" = S) and completely interrupting the passage 
of supercurrents in this wing. The passage of any nor- 
mal component of current density is made possible by 
filling up the cut with a normal metal N . The carrier 
type of the normal-metal insert N was chosen to be ei- 
ther an electronlike N^ (silver) or holelike Nh (indium). 
Thermal voltages resulting from the same temperature 
gradient acting on both continuous and normal-metal- 
filled halves of the sample were detected as a function 
of temperature around Tc. The measured difference be- 
tween the thermopowers of the two halves AS* = Sr — Sl 
was found to approximately follow the linear dependence 



A5(T)~5p-S(Tc-T), 



(1) 



where Sp = AS'(Tc) is the peak value of AS'(T) at 
T = Tc, and i? is a constant. The best fit of the ex- 
perimental data with the above equation yields the fol- 
lowing values for silver (Ag) and indium (In) inserts, re- 
spectively: (i) Sp{Ag) = ~0.26±0.01fiV/K, B{Ag) = 
-0.16 ± O.lfiV/K'^; (ii) Sp{In) = 0.83 ± O.OlfiV/K, 
B{In) = Q.n ±Q.l^Jy/K^. 
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FIG. 1. Schematic view of the sample geometry with 
S''Ai'5'"-junction and contacts configuration. The thermopow- 
ers Sr and Sl result from the thermal voltages detected by 
the contact pairs 4 — 5 and 1 — 7, respectively. 



B. Interpretation 

Assuming that the net result of the normal-metal in- 
sert is to break up Cooper pairs that flow toward the 
hotter end of the sample and to produce holelike (In) or 
electronlike (Ag) quasiparticles, we can write the differ- 
ence in the generalized GL free energy functional /S.Q of 
the right and left halves of the C-shaped sample as 
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(2) 



2. N — Ag (electronlike metal insert) 



where 



ATiij] ^TR-TL = a{TM' + ^\^l;\^ (3) 



and 



(4) 



Here tp = j-^le"^ is the superconducting order parameter, 
fip and fiq are the chemical potentials of quasiparticles 
and Cooper pairs, respectively; a(T) = a{T — T^) and 
the GL parameters a and (3 are related to the critical 
temperature T^, zero-temperature gap Aq = l.TB/csTc, 
the Fermi energy Ep^ and the total particle number den- 
sity n as a = 2AokB/Ep and (3 = aTc/n. 

As usual, the equilibrium state of such a system is de- 
termined from the minimum energy condition dQ/d\^\ = 
which yields for T < 



(5) 



Substituting IV'oP into Eq.(2) we obtain for the general- 
ized free energy density 

Ar!(r) = AGi^jo] : 



[a{T, -T) + A^i]^ 



2(3 



(6) 



In turn, the observed difference of thcrmopowers AS(T) 
can be related to the corresponding difference of trans- 
port entropies Act = -DAfl/dT as AS{T) = Aa{T)/nq, 
where q is the charge of the quasiparticle. Thus finally 
the thermopower associated with a pair-breaking event 
reads 



ASiT) 



Afj. 
qTc 



2AokE 



(7) 



where Ep = Ep — fiq accounts for the shift of the Fermi 
energy Ep due to the quasiparticle chemical potential 
/ig . Let us discuss now separately the case of In and Ag 
normal-metal inserts. 



1. N = In (holelike metal insert) 

In this case, the principal carriers are holes, there- 
fore g = -fe in Eq.(7). Let the holelike quasiparticle 
chemical potential (measured relative to the Fermi level 
of the free-hole gas) be positive, then fiq ~ +fi and 
A/i = fjLq — fjip = fi — (— 2/i) — 3/i. Here /ip — — 2/i 
comes from the change of the pair chemical potential of 
the holelike condensate with respect to the holelike quasi- 
particle branch. Therefore, for this case Eq.(7) takes the 
form 



AS\T) = i(^^' ^ 



knTr 



2Aofc£ 



{Tc-n (8) 



The principal carriers in this case are electrons, there- 
fore q = —e. The electronlike quasiparticle chemical po- 
tential (measured relative to the Fermi level of the free- 
hole gas) is — /i. Then jiq = — /i and Aji = —ji — {—2fi) = 
fi. For this case Eq.(7) takes the form 



AS'^(T) = - — 



fc_B 

e 



2AofcB 
eE^pT, 



(Tc-T), (9) 



where Ep = Ep + fi. 

Using the above-mentioned experimental findings for 
the slope B and the peak Sp values for the two normal- 
metal inserts (see Eq.(l)), we can estimate the order 
of magnitude of the Fermi energy Ep and quasipar- 
ticle potential /x. The result is: Ep = QAQeV and 
fi = 5 X 10~^eV, in reasonable agreement with the other 
known estimates of these parameters. Besides, as it fol- 
lows from Eqs.(8) and (9), the calculated ratio for peaks 
\Sp{In)/Sp{Ag)\ = 3 is very close to the corresponding 
experimental value IS^'Ziln)/ S^''P{Ag)\ = 3.2 ± 0.2 ob- 
served by Gridin et al llSl . 



III. SND CONFIGURATION: PREDICTION 

Since Eqs.(2)-(4) do not depend on the phase of the 
order parameter, they will preserve their form for a 
DND junction (created by two d-wave superconductors, 
S' = S" = D, see Fig.l) bringing about the result similar 
to that given by Eqs.(7)-(9). It means that the experi- 
mental method under discussion (and its interpretation) 
can not be used to tell the difference between SNS and 
DND configurations, at least for temperatures close to 
Tc- As for low enough temperatures, the situation may 
change drastically due to a markedly different behavior 
of s-wave and d-wave order parameters at T <C Tc- As 
we will show, this method, however, is quite sensitive 
to the mixed SND configuration (when S' = S has an 
s-wave symmetry while S" = Z? is of a d-wave symme- 
try type, see Fig.l) and predicts a rather specific relative 
phase {9 = (f>s — (f>d) dependences of both the slope B{0) 
and peak Sp{9) of the observable thermopower difference 
AS{T,e). 

Following Feder et al who incorporated chemical 
potential effects near twin boundaries into the approach 
suggested by Sigrist et al , we can represent the gen- 
eralized GL free energy functional AQ for SND configu- 
ration of the C-shaped sample in the following form 



Ag[^Ps,H = ^Gii^s] + AGli^d] + Ag,, 



where 



AGiibs] = AT[ibs] - AVIV'S 



(10) 



(11) 



where Ep = Ep — /i. 



AGl^Pd] ^ ATlijd] - Afi\^d\\ 



(12) 
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and 



AS{T,e) = Sp{9)-B{9){T,-T), 



AG,, 



Tii^.nv'dp + f (v-rv^+^v-f) (13) 



Here ipn = is the n-wave order parameter 

(n = {s, c?}); AjF[-!/)s is given by Eq.(3) with the corre- 
sponding parameters as{T) — as{T — Tcs),l3s, ad{T) — 
ad{T — Ted), and f3d for s-wave and d-wave symmetry, 
respectively. 

An equihbrium state of such a mixed system is deter- 
mined from the minimum energy conditions dQ/dltpsl ~ 
and dQ/d\ipd\ — which resuh in the foUowing system 
of equations for the two equihbrium order parameters 
and tpdo 



^(T)IVsol 
Ad{T)\^do\ 



-/3d IV', 



do 



-mi^PdoiM 



where A,{T) = a,(r) - A^, AdiT) 
we introduced relative phase = (ps 
rameters 



m = 71 - 

A{9) = 5i 



72 cos 29 
- 82 cos 9 



A(0)IV'do| (14) 
= A(^)|^.o| (15) 

ad{T) — Afj,, and 
(pd dependent pa- 

(16) 



Notice that the A{9) term favors 9 = In {I integer), while 
the T{9) term favors 9 = nn/2 {n — 1, 3, 5 . . .) which cor- 
responds to a T-violating phase p3[ |. In principle, we 
can resolve the above system (given by Eqs.(14)-(16)) 
and find '0„o for arbitrary set of parameters a„, /3„, and 
Ten- For simplicity, in what follows we restrict our con- 
sideration to the two limiting cases which are of the most 
importance for potential applications. 



A. Twin boundaries in orthorhombic d-wave 
superconductors 

First, let us consider the case of similar superconduc- 
tors comprising the SND junction with IV'sol = IV'dol = 
IV'ol, as = ad = a, Ps = Pd = /?, and Tcs = T^d = T^. 
This situation is realized, for example, in a d-wave or- 
thorhombic sample (like YBCO) with twin boundaries 
which are represented by tetragonal regions of variable 
width, with a reduced chemical potential |]l3| . In this 
particular case, Eqs.(14) and (15) yield for T <Tc 



a{Tc ~T) + An + A(6>) 
f3 + T{9) 



(17) 



After substituting the thus found jV'ol into Eq.(lO) we 
obtain for the generalized equilibrium free energy density 



An{T,9)=Ag[ijo] 



[ajT, -T) + A^l + A{9)]^ 

p + r{9) 



(18) 



which in turn results in the following expression for the 
thermopower difference in a C-shaped sample with SND 
junction (see Fig.l) 



where 



and 



SpiO) = - 



2Afi / 1 + ^cos( 



qTc V 1 + 7 cos 29 



B{9) = 



4AofcB 
qEpTc Vl+7cos20 



(19) 



(20) 



(21) 



Here, 7 = 72/(71 + /?) and S = 82/ {5i + Afi) with A/i 
and P defined earlier. 




FIG. 2. Predicted phase-dependent thermopower response 
of SND configuration in a C-shaped sample (see Fig.l). Solid 
and dashed lines depict, respectively, the relative phase 9 de- 
pendence of the normalized slope B{9) / B{Q) and peak value 
Sp{9) / Sp(0) of the induced thermopower difference, according 
to Eqs.(20) and (21) with 7 = 5 = 1/2. 

Fig. 2 shows the predicted 6'-dependent behavior of the 
normalized slope B{9)/B{0) (solid line) and the peak 
Sp{9)/Sp{0) (dashed hne) of the ^A^ZJ-induced ther- 
mopower difference AS{T,9) near Tc, for •y — 6 — 1/2. 
As is seen, both the slope and the peak exhibit a max- 
imum for the s ± id state (at 9 — 7r/2) and a minimum 
for the s — d state (at = tt). Such sharp dependen- 
cies suggest quite an optimistic possibility to observe the 
above-predicted behavior of the induced thermopower. 
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using the described in Section II sample geometry and 
experimental technique. Besides, by a controllable varia- 
tion of the carrier type of the normal- metal insert, we can 
get a more detailed information about the mixed states 
and use it to estimate the phenomenological parameters 
7i,2 and 5i^2- 



B. Low- Tc s-wave superconductor and high-Tc d-wave 
superconductor 

Let us turn now to another limiting case when the two 
superconductors of the SND junction are different, so 
that IV'sol ^ IV'dol, a-s ^ ad, (is ^ Pd, and T^s ^ Ted 
but the charge imbalance effects are rather small and 
can be safely neglected, that is we assume A/x = 0, and 
A(9) = 0. Such a situation can be realized for an s-wave 
low-Tc superconductor (like Pb) and a d-wave high-Tj, 
superconductor (like orthorhombic YBCO) In fact, 
the solution for this particular case is well-known. It 
has been discussed by Sigrist et al in a somewhat 
different context. The corresponding expressions for the 
equilibrium order parameters read 



1^, 



sO\ 



'do I 



pdasiT)-Tie)adiT) 
T^9)~l3s(3d ' 

/3,arf(r)-r(g)a,(r) 

T^{9)-(3s/3d ' 



(22) 



(23) 



where as{T) = as{T - T^s) and ad{T) = ad{T - T^d)- 

After substituting this solution into Eq.(lO) we obtain 
for the the thermopower difference 

A5(T, 9, t,) = Spi9, t,) - B{9, U){T,s - T), (24) 

where tc = Tcd/Tcs and both the peak 

Spi9,Q^{l-tc)f{9,te) (25) 

and the slope 



B{9,tc) = 2- 



are governed by a universal function 



qEF[ri9) + Pstc 



(26) 



(27) 



Notice that in this case (when changes in chemical poten- 
tials can be neglected) the peak's amplitude Sp{9, tc) will 
be entirely dominated by the critical temperatures differ- 
ence Tcd — Tcs of the two superconductors because f{9, tc) 
is a smooth function of tc- It would be interesting to test 
the predicted behavior of the induced thermopower in a 
C-shaped sample with an SND junction (see Fig.l) using 
a low- Tc s-wave and a high- Tc rf-wave superconductors. 



In summary, to probe into the mixed s±id pairing state 
of high- Tc superconductors, we calculated the differen- 
tial thermopower AS of SND junction in the presence 
of the strong charge imbalance effects (due to a nonzero 
difference between the quasiparticle fiq and Cooper pair 
/ip chemical potentials) using the generalized Ginzburg- 
Landau theory for a homogeneous admixture of s-wave 
and d-wave superconductors near Tc. The calculated 
thermopower was found to strongly depend on the rel- 
ative phase 9 = (ps — 4'd between the two superconduc- 
tors exhibiting a pronounced maximum near the mixed 
s i id state with 9 — ±7r/2. The experimental condi- 
tions under which the predicted behavior of the induced 
thermopower could be observed were discussed. 
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